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'■ Abstract 

qh| The so(2, 1) Lie algebra is applied to three classes of two- and three- 

dimensional Smorodinsky-Winternitz super-integrable potentials for which 
the path integral discussion has been recently presented in the literature. We 
have constructed the Green's functions for two important super-integrable 
5^ I potentials in R^. Among the super-integrable potentials in R^, we have con- 

sidered two examples, one is maximally super-integrable and another one 
minimally super-integrable. The discussion is made in various coordinate 
systems. The energy spectrum and the suitably normalized wave functions 
of bound and continuous states are then deduced. 
PACS 03.65-Quantum theory ; quantum mechanics 
PACS 03.65.Fd -Algebraic methods 
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I. Introduction 



Ever since the success of the algebraic approach based on the non-compact 
groups in calculating the wave functions and the transition amplitudes for 
the hydrogen atom by Klcincrt^ , a renewed interest for this method has 
been emerging. Hence a large amount of physical problems for which the 
Schrodinger equation reduces to the confluent hypergeometric equation have 
been treated in the framework of this approach. In particular, the radial 
harmonic oscillator, the radial Coulomb and Morse oscillator potentials and 
the Natanzon confluent potentials^ , which generalize the latter, have been 
discussed in a variety of SO (2, 1) realizations^"^ . 

The above so(2, 1) algebraic approach has also been applied to noncentral 
potentials such as the ring-shaped potentiaP° introduced by Hartmann^^ to 
describe the molecular interaction of cyclic polyenes and the so-called double 
ring-shaped oscillator^^ which is the Quesne ring-shaped oscillator^^ plus an 
r^^ sec^ 9 term. 

The (MS) variant of the algebraic method summarized in section II has 
been developed by Milshtein and Strakhovenko^^ to construct the Green's 
function associated with the problem of Dirac electron in a static Coulomb 
field. The straightforward algebraic calculation of the Green's function for a 
given potential represents an undeniable advantage which allowed this variant 
to be given a great deal of attention, in recent years, following the develop- 
ment of path integration techniques. Hence a remarkable set of potentials 
has been studied in this algebraic approach. Among them, we can quote 
potentials of practical interest, such as Morse's^^ , the radial harmonic os- 
cillator and the radial Coulomb potentials . The Hartmann ring-shaped 
potential^^ , the compound Coulomb plus Aharonov-Bohm potential^'' and 
two highly singular non-isotropic potentials associated to a highly distorted 
spherical Coulomb field with an additional double ring well and a highly 
distorted cylindrical Coulomb field have also been studied in parabolic co- 
ordinate systcms^*^ . Within the framework of the to non-bijective 
Kustaanheimo-Stiefel mapping, the Kaluza-Klein monopole system^^ and a 
noncentral potentiaP° which generalizes the Coulomb potential and the Hart- 
mann ring-shaped potential and also, due to its close link with the latter, 
the compound Coulomb plus Aharonov-Bohm potential have been treated 
by means of the same algebraic approach as Milshtein and Strakhovenko. 
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Cylindrical parabolic coordinates have also been used in discussion for the 
so(2, 1) algebraic method of another type of noncentral potentials^^'^^ . All 
these potentials discussed with the help of the (MS) variant appear in the 
general classification of potentials in two and three dimensions possessing 
dynamical invariance groups initiated about 30 years ago by Smorodinsky 
and co-workcrs^'^ , continued by Kiblcr and Winternitz^'^ , and revived, in 
recent years, by Evans'^'^''^^ . This classification was established according 
to the number of degrees of freedom, quadratic integrals of motion in the 
momenta and coordinate systems in which the potential allows the separa- 
tion of variables. The Hamiltonian systems with these potentials are called 
super- intcgrable. Generally, in n dimensions, a system is called "minimally" 
super-integrable if it has (2n — 2) constants or integrals of motion ( includ- 
ing energy ), and it is called "maximally" super-integrable if it has (2n — 1) 
integrals of motion^^ . A list of minimally super-integrable and maximally 
super-integrable potentials with the corresponding constants of motion in the 
classical form and all separating coordinate systems has been established by 
Evans^^ . On the basis of this classification, Grosche et aP^ have recently 
presented a detailed path integral discussion of the so-called Smorodinsky- 
Winternitz super-integrable potentials in many coordinate systems. It is 
to be noted that almost all the potentials contained in this classification in- 
volve centrifugal or angular barriers which possess point singularities. Conse- 
quently, following Kleinert^'' the time sliced path integral for these potentials 
does not exist in any coordinate system. So, it is necessary to regularize the 
system in question by an appropriate set of new coordinates in order to find a 
path integral expression without collapse. This problem does not occur in the 
framework of the (MS) variant of the so(2, 1) algebraic approach owing to its 
local (differential) character. This constitutes a great advantage in studying 
the problem of the Smorodinsky- Winternitz potentials in this framework. 

The plan of this article is as follows. We briefly review the so (2, 1) Lie 
algebra and its use in calculating Green's functions in section II. We study 
a set of two potentials in two dimensions, algebraically obtain the Green's 
function in the various separating coordinate systems and deduce the en- 
ergy spectrum as well as the corresponding normalized wave functions in 
section III. Sections IV and V deal with two examples of three-dimensional 
potentials. The construction of the Green's functions is made in different co- 
ordinate systems. The energy spectrum and the normalized wave functions 
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are evaluated. Section VI will be a conclusion. 



II. Green's function and so(2, 1) Lie algebra 

Let's briefly review the main features of so(2, 1) Lie algebra and its use 
in the calculation of the Green's functions to make this paper self-contained. 
A set of three operators {Ti,T2,T3} characterized by the commutation rela- 
tions^^ define it: 

[Ti, T2] = -zT,, [T2, T3] = -tTs, [Ti, T3] = -tT^. (1) 

Because of the type of potentials we shall deal with in this paper, we have 
to use the following differential realization of the operators: 



^2 fd^ 



d 1 



M 



2MW r^^^"^^"n"9^^2j'^^^"^ = 4n^""' 

(2) 

with < a; < 00. 

By using Schwinger's integral representation'^^ , the Green's function as- 
sociated to a potential V{x) with the S0(2, 1) group symmetry is given by 



/■oo 

G{x,x']E)— / dS exp 
Jo 



K{x, x'; S), 



(3) 



where 



K{x, x'; S) — exp < — 



tS 
iS 



2M 



+ V{x) 



6{x - x') 



exp <! — ^ [Ti{x) + 2hWTs{x)] \ S{x - x'). (4) 



The calculation of this kernel is based upon the use of two Baker-Campbell- 
Hausdorff formulas^^ 



exp I [Ti + 2hWT3\ \ = exp(-iar3) exp(-i6r2) exp(-icri), (5) 
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where 



and 



with 



a — 2hu}ta,n(u}S), b — 2Ln [cos(a;5')] , c — tan(a;5'), 

nw 



exp(— iaTa) exp(— i/3T2) exp(— i7Ti) = exp(— icTi) exp(TT3) 



(6) 



(7) 



a 



1 - iTc/2 



, p^2Ln[l- 



ITC 



7 



l-iTc/2' 



(8) 



Here, we also have to use the Laplace transform of the Dirac distribution 



5{x - x') 



2n? 2m 



ioo+5 



drexp 



ioD+5 



; <^<o, (9) 



in order to obtain a manageable result as follows: 



x" = x". 



exp(-ijTi)x^ = 

Using relations (9), (5) and (7), the kernel (4) can now be written 



(10) 



K{x, X \ S) 



M 

2h? 



x'^ exp ( — mTa ) exp {—ibT2 ) x'' cxp 



[ 2h?c 



rzoo+S exp(-^Tlexp(;^^^l 

X / dr ^ ^ 



where the well-known formula 

exp(-i/3r2)/(x) = exp 



(11) 



(12) 



has also been used. 

The integral can be calculated thanks to the residue theorem after the 
exp ^^p^ ^_)_2j/c ) series has been effected. Hence, we obtain 



K{x,x']S) 



ihsm{ijjS) 
^ , Mojxx' 



XX' exp 



[X +x )cot(c<;6) 



ihsm{u!S) J ' 



(13) 
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where A = — | and I\{x) is the modified Bessel function. 

We can now use (13) for any coordinate in a multi-dimensional system, 
provided that the {H — E)~^ inverse resolvent operator can be transformed 
into a linear combination of the above mentioned Ti{i — 1, 2, 3) operators. 

III. Two-dimensional maximally super-integrable po- 
tentials 

We shall study here a set of two important potentials belonging to a class 
of two-dimensional Smorodinsky-Winternitz potentials. They are character- 
ized by the existence of three functionally independent integrals of motion, 
which means that there is a pair of quadratic operators corresponding to 
these integrals of motion which commute with the system's Hamiltonian. 
The number of such integrals being superior to that of degrees of liberty, 
they are thus called maximally super-integrable potentials. 



A. Let us study the potential 



^/xl+xl 4M ^yxf + xl \y^xl+xl + Xi y^xf+xl — Xi 

with positive Oq, ki and k2 constants . It admits the following three func- 
tionally independent integrals of motion: 



T ^ {L3,P2} , a(r?-0 _ 2(fc^-l) 2(fc2-l) 
-'2 4M 



■-1 



(15) 



This potential is exactly solvable in two coordinate systems, namely 
parabolic and polar. For ki = k2 = \, equation (14) can be reduced to 
the Coulomb potential treated with a path integral approach"^"^'^^ . The al- 
gebraic solution to this potential via the (MS) variant is easier to establish 
using the parabolic coordinates or even the Levi-Cevita variables^^ defined 
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by xi = uf — U2, X2 = 2uiU2 (—00 < ui,U2 < +00). The Green's func- 
tion G(7?, p' E) associated to the potential (14), in Schwinger's integral 
representation, is given by 



3 4 



Ap^ \ 2M law? 

2 



(16) 



By applying the time transformation 5* — > r defined by r = ^, the 
Green's function (16) can be written 



G{ p , p' ; E) = I drexp 
'0 



-{4ao + iO)T 



1[K{u„u'.,t), (17) 



where 



K{uj, u'-, t) = exp \ -J [T,{uj) + 2fi'u;''n{uj)] \ 6{u, - u'^), (18) 



8E 
' M ■ 



with /ij = I ± kj and a; = 

Using the equation (13), we can write the propagators (18) as follows: 



Muj ^Uju'- 
K{uj, M ■; r) = . — - exp 
zAism(a;r) 



'iMuj 



2h 



MujUju'j 



(19) 



1. Parabolic coordinates 



In the parabolic coordinates ^ = \/2ui, r] = \/2u2 and after changing 
a; — > 2a; and t | and by taking into account (19), the Green's function 
(17) becomes 



G{~p,p']E)— I drexp 
Jo 



-{2ao + i0)r 

f b 



m,V,e,v';r), (20) 
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where 



Muj 



ihsm{uJT 
iMuj 



with 



X exp 



2h 



ie + V' + e + V") cot(a;r) 

Mujrjri' 



i^sin(a;T) / ' 



2E 



(21) 



(22) 



To find the energy spectrum and the normahzed wave functions of the 
bound states, we make use the Hille and Hardy formula^^ 



1-z 



—z- 



1-z 



z \< 1, 
(23) 



where the Lf^{x) are the Laguerre polynomials, and the integration over S 
yields the quantization condition 

ni+n2±y ±y +p = 0. (24) 
where p = — Therefore, the Green's function (20) can be written as: 



,712=0 



E + iO- En„n, 

with the normalized wave functions (a — is the Bohr radius) 

2 ni!n2! 



(25) 



^ni,n2 (C) V) 



o?N^ r(ni ± A;i + l)r(n2 ± ^2 + 1) 

X±h2. 



aN 



2 X 



oN 



exp 



2aN 



e 



aN I \aN )' 

(26) 
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and the discrete energy spectrum given by 



^;»..« = -^i /V = n, + n,±|±| + l. (27) 

To determine the energy spectrum and the wave functions of the continu- 
ous states, let's go back to the expression (21) and use the dispersion formula 
(see Ref.37 , p. 884, Eq. (7.694)) 



exp [-(x + I/) cot aj /2;^ 



smo; V sma 



Given x = ff V = and a = a;T , we obtain 
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^fl ki \^ fl ki 

xr ± I + zp,) r Q ± I - ip, ) e(---^)(^^+^^) 

x*^-.p„±i M.,,,,^ (^f 

x^-.,,,** (^"^j A'.,„.^ (^-J'^j . (29) 

If we now transfer (29) into (20), we shall obtain the poles of the continuous 
state Green's function by integration on the r variable. They will be defined 
by 

^(p«+P.)-^ = 0. (30) 
If we convert this into energy via (22), the values of the energy will be 

-ilil with p=—^ r. (31) 
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Let's now change the variables defined by 



In this case, we can write the Green's function as follows: 



POO roo 

G%C, V, e, v'; E)^ih dp d<; 

Jo J-oo 

with the wave functions given by 



E + iO-Ep 



2p\a^^)'^ 2 



27rr(l±A;i)r(l±A;2)v^ 

fc2 



d-W^^)M_^,._.^^{-tpV^). 

2p\a ^''^ 2 



2. Polar coordinates 

Let's use the (p, 0) polar coordinates defined by 



^ = \/2mi = \/2pcos — , 77 = V2u2 

and the addition theorem formula^^: 

-/i,(2;sinQ;sin P) I /^{z cos a cos P) 



2/9 sin — , 
^ 2' 



(32) 



(33) 



g 2ap 



(34) 



(35) 



= (sin Q;sin/3)''(cosa cos/3)'^ > (i/ + + 2n + 1)- 

r(i/ + n + l) 

n=0 ^ ' 

r(i/ + /i + n + 1) 



X- 



r(// + n + 1) 



^W2n+i(^)P^'^Hcos(2a))Pi^''^)(cos(2/3)). 



(36) 



This will give us the following form of the Green's function (20): 

G(y , 7; i^) = E Gr.{p. p; E)$(^'=-^^^) (0 ^i^'^-^^^) (f ) , (37) 
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where the angular wave functions are those defined in function of the 
p^±fc2,±fei)^^^^ jg^gQ|-,j polynomials by 



X (sin 0) si'^^ (cos (f))^^''' P^^'^^.ifci) (cos(20) ) . 
The radial Green's function Gn{p,p'', E) included in (37) is defined by 



(38) 



Mlj 



dr 



ih Jq sin(ci;r) 
'iM(jj 



exp 



-{2ao + i0)r 

I V 



X exp 



2h 



{p + p') COi{ijJT) 



^ f2Mu^fpi7 
\i^sin(a;r) 



(39) 



with A = n + 1(1 ± fci ± k2). 

Then, thanks to formula (see Ref.^^ , p. 729, Eq. (6.699.4)) 



oo g-2pg 

smh q 



r(27 + i)v^ 



exp — [x + y) coth q 



'27 



sinh q 



(40) 



valid for Re(p + 7 + |) > 0,Re(7) > and y > x, where M_p^^(a;) and 
W-p^^{y) are the Whittaker functions, we can write (39) as follows: 



2iu!T[2X + IjvPP 



/2MCC) ; 



/2Mcu 



-p 



(41) 



and p > p'. 



Prom (41) and (37), we can deduce that the complete Green's function is 
given by 



/2Muj \ 



r(p + A + |) 

2iuJ^{2X+l)^/^ 



(42) 
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The normalized wave functions and the energy spectrum of the bound 
states are given by expression (23) apphed to Gn{p, p'; E) as defined in (39) 
provided we make the adequate change of variables and thus obtain: 



*m,n(p, 0) = 



ml 



2p 



a{m + \ + \) 



a(m + A + I) / \a{m + \+\ 

t 
2^ 

In order to evaluate the contribution of the continuous spectrum to the 
Green's function, let's write (41) as follows: 

ih f dz , 1, 



A7rT{2X + l)y/pp' Jc E + tO- ^ 2 



where C is the closed contour, 

r z = k;ke [-R, R] . . 

^ = i?e*^0e(7r,27r). ^^^^ 

At the i? — > oo limit, taking the asymptotic behaviour of the Whittaker 
functions (see Ref.^''' , p. 1061, Eq. (9.227)) into account, it is easy to 
demonstrate that the integral over the semicircle vanishes. Which leads to: 

= 4.r(2A + i)V^y_£; + io-|r'^^ + '+2) 

xM_p,;,(-2iV)W^-p,A(-2iA;p)- (47) 

By using the following formulas (see Ref.^^ , pp. 1061-2, Eqs. (9.231.2) and 
(9.233.1)) 

M^,^{z) = e-'<^+"^'^M_^,^{-z), with 2p ^ -1, -2, -3, (48) 
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and 



r(/x - A + 



(49) 



valid for aigz e ] — | [, and 2// 7^ —1, —2, —3, the expression (45) can 
be written 



ih 



47rr2(2A+l)Vpp^ Jo 



dEk- 



E + iO-Ek ' 



(50) 



where Ek 



2M 



and the radial wave functions given by 



M 



|r(p + A 



1 \ I HP 

n) e 2 



r(2A + i) VP 



(51) 



B. Let us now study the potential 

Let us now study the potential 

^xl+xl+xi+fi2\l \/x{^. 



V2{P)- 



+ 



l-Xi 



(52) 



b2 Y ^i-r^2 

with real j3i and constants. This potential has the following three func- 
tionally independent integrals of motion: 

TT _ P2 T , J _ {L3,Pi} a(g-r?)+/3ir?.yg72-/32$VV2 

-"2 - 2M + *^2l P j, -'1 - 1+^ /roN 

_ {Lz,Pi} _ (^jy/i/i-gy^) ^ 

For given parabohc coordinates Xi ~ |(^2 — rf), X2 = ^t], > and 
7] G R), in the Schwinger's integral representation, the Green's function 
associated to this potential can be expressed 



G{p,p';E 



dS exp 
X exp 



■f(i. + .o) 



2p 



(54) 
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where 



2M2p \de dr]y p p 



with p = U^^ + v'^). 



If we now perform the time transformation defined by a = ^ = 



LDC 



to separate the ^ and i] variables and use the mutually orthogonal parabolic 
coordinates (^^^V^ ~^ ~ ~ arrive at 



g{p,p']E]= dcrexp 



E 



(56) 

where each of kernels K{^, a) and K{ri, rf; a) can be treated with the 
so (2, 1) Lie algebra and we have 

K{x,x';a) = ^ J] exp |-|^ [ri(x) + 2;iVr3(x)] 1 5(x - x') 



//=0 



Muj 



■exp 



iMuo 

2h 



X +x ) cot{u;a)- 



2xx' 



sin(a; cr)_ 



(57) 



with X = or 77^ and cu — 

Substituting (57) into (56), we obtain 



2E 
' M ■ 



G[p,p';E] = — - / ■ , , exp 



xexp 



2inh Jo sm(uja) 

{iMuj 



J [2a. 



13! + PI 



E 



sm(uja) 



(58) 



In order to determine the energy spectrum and the normalized wave func- 
tions of the bound states of the physical system, let's apply the Mehler for- 
mula^^ 



1 



: CXp 



1 



2(1 -a2) 



[{x"^ + x'^){l + a^) - Axx'a] 



14 



exp 



oo 
n=0 



(59) 



With the help of an adequate change of variables, the poles of the Green's 
function (58) will be obtained thanks to an integration over a; the discrete 
energy spectrum is found by solving the equation 



OJ^ - -TT^uj^ - 2 ' l_ = 0, or as well En ^ 



Nh 



NMh 



(60) 



with iV = ni + 712 + 1. CKo being positive, then this cubic equation has one 
real root: 



UN 



2qo 
3Nh 



+ Ai + As, 



where 



A, 



2aoV , f^f+f^l 



3NhJ NMh 



NMh \ NMh 



+2 



2ao 
3Nh 



(61) 



(62) 



with {j = 1,2). 

We may obtain the normalized wave functions of the bound states from 
the residues of the integrated expression of the Green's function (58), 



M 



ni,n2 



1 



\ ^ ^ ^ Nh^ ^ NMh 



lim 



xexp 



Muj 



2h 



H„ 



Mujn 



h 



Muon. 



h 



(63) 



Here, it is obvious that only states with an even total number of oscillator 
quanta contribute. 

To find the wave functions of the continuous states, let's go back to the 
Green's function (58) and make use of the following relation 



^/2 



TTsmo; 



: exp {—{x-\-y) cot a) exp 



2^ 



sm a 



15 



(27r)^ 



X 



+ 



^^g{,r-2a)p 



+ ip 



(64) 



which is estabhshed from the dispersion formula (28). The Ei^\z) and 
are even and odd parabohc cyhnder functions with respect to the 
variable z, respectively^^ . The poles of the continuous state Green's func- 
tion will be obtained by integration on the r variable. They will be given 

by 



+ Pi 
E 



0. 



Then, by performing the change of variables 



1 A 



2p 



where a 



M{ao-M[0l+Pl)/t9p^) 



(65) 

, it is possible to 



write the contribution of the continuous part to the Green's function as 

G\i. ri. e, V'; E) = ih r dp r (66) 



where Ep — 
form 



and the continuous functions '^p,<;{^-,r]) have the 



V r (t+i Ci+^)) (-^-'-'wl 



\,r(f+^g-.)jBi'i^,(._^)(e-/4V2pj) 



(67) 
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IV. Three-dimensional maximally super-integrable po- 
tentials 



In three-dimensional Euclidean space, Smorodinsky and co-workers have 
found a set of five potentials which have five functionally independent in- 
tegrals of motion. These three-dimensional potentials are called maximally 
super-integrable potentials. At least, each potential of the so-called class of 
Smorodinsky- Winternitz potentials can be treated in two coordinate systems 
through the so(2, 1) Lie algebraic approach. Here, we shall restrict ourselves 
to study the potential 



r = h 

' r 2M 



+ 



XT 



The integrals of motion are 



h-- 



- jl- 

2M 
2M 



2M 



+ 



2M tan^ 9 cos^ 4> ' 



2M \ cos^ (f, 
2M 



X2 



sin <p 



(68) 



T _ 1 



2Msin^6» 

n? 



cos^ 
2 



sin*^ <t> 



(69) 



1x1X2 + Ix2Xi) (^ V) 2M^?? (cos2^ + ih?^) 



where 1^ = {Lj, Pj} = LiPj + PjLi with = {xi, X2, X3). 

This potential is a generalization of the Coulomb potential analyzed by 
various authors in the path integral*^^"^^ and algebraic approach^^'^^ contexts. 
As it features singularities for xi = and X2 = 0, all we need is to discuss it 
in the < xi, a;2 < 00 and X3 E R area. It is possible to evaluate the Green's 
function 



GiT", r';E) 



dS exp 



E-iO) 



5(T* — r' ) 



(70) 



in parabolic and spherical coordinates. 



1. Parabolic coordinates 



For given parabolic coordinates Xi — ^77 cos 0, X2 — ^77 sin 0, 
^3 = |(C^ ~ ^^)) ^ > ^iid < < 27r, the Green's function (70) can be 
written 



^ /-oo 

G{'r',r']Ey dS exp 
Jo 



_!|(/J3_£;_i0) 



, (71) 
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with 



2M 
2aQ 



Id Id 



1 



^2 _|_ ^2 ^(5^2 ^ Qj^2 Tj drj J ^"^Tf d(i)^ 

^2 /U2 1 u2 1 



+ 



2 1 V.2 

4 , '^2 4 



(72) 



^2 _|_ jy2 2M ^^7^2 \^ (.Qg2 gjj^i 

Separating the angular part of the expression (71) by a time transformation 
s 



defined by r = we can deduce that 



G(7^, r',E) 



(73) 



where 



i^(e,^,r,V;T) = exp<|-|[i73(e,^)- W-iO] 



X 



(74) 



and 



K{(t),(t)'-T) = exp 



sm 



5(0-00 



^ 2M V cos2 , 

$p.±^O(^)$(±^.±^O(^0 exp f-^||r ) , (75) 



n=0 



n 2M 



with Ai = 2n ± /ci ± A;2 + 1- 

It is to be noted that the exphcit construction of the kernel K{(}), 0'; r) 
thanks the (MS) variant of the algebraic technique is being investigated and 
will be the subject of our forthcoming publication. 

Introducing (75) into (73) and applying the r ^ S inverse time transfor- 
mation will allow us to write (73) as follows: 



^ °° poo 

G{-T',r';E) = ^$p2'±'=i)(0)$i±'=2'±'=i)(0') / dSexp 
n=o 



'iS_ 



(E-iO) 



iS 

xexp^-- 



] sj^ - - V) 



(76) 
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where 



^2 Id Id 



2ao tfXl 

If we now eliminate the |^ and operators by applying ^ "I" |^ ^'^'-^ 
and "^^J^-* respectively, we can then proceed with a new 
S ^ a time transformation defined by cr = |2^^, which will allow us to write 



drp Tj dr) 



n=0 



with 



:exp^-- 



2M \du^ 



A? 



2 1 
4 



-Eu' 



-{2ao+iO)a 



(78) 

(79) 



where w = (^ or 77). 

So, it is possible to give the kernel (79) in function of the so(2, 1) Lie 
algebra operators. Indeed, by following the equations (2) and (13), we obtain 



1 



h 

'iMuj 



Kniu, u'; a) = exp <^ -— \Ti{u) + 2hWn{u)] } 5{u - u') 

Vuu' \ n ^ -I I 

Mu 



ihsiniuja 



-exp 



2h 



u +u' )cot{u;a] 



Mujuu' 



ihsiniua 



with uj = sj—^- The Green's function (78) can now be written 

00 

G(r,7;E) = ^$(±'=-±^^)(0)$(±'=-±^i)(</.')G„(e,r/,e',r;';E), 

with 



(80) 



(81) 



n=0 



G„(e,r^,e',V;E) 



Muo 
ih 

X exp 



da 



sin (coa) 



exp 



h 



{2ao + iO)(T 



2h 

ihsmiuja] 



Mujrjri' 



(82) 
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Thanks to the Hille and Hardy formula (23) for the discrete part and 
thanks to the scattering relation (28) for the continuous part, the Green's 
function (81) can be developed into partial waves as follows: 



G{r ,r = ^^2^< 2^ 



n=0 \ni,n2=0 
poo poo 

/ 

Jo J- 



E + M-En 
*,,«,„(e,r^,0)vi>* (r,V,</'') 



+ , - I .(83) 

Hence, for bound states-the normalized wave functions and energy spectrum 
will be 



ni\n2\ 



a?N^V{ni+ Ai + l)r(n2+A2 + l) 



2 

aN 



X exp 



2aN 



ni 



aN 



(84) 



E. 



Ma, 



""^'Wn^' Ar = ni+n2 + Ai + l, (85) 

and, for the continuous states, the normalized wave functions and the energy 
spectrum will be 



27rr2(Ai + l)irj 



g 2ap 



^(±fc2,±fcl)j-0-J^ 



2M 



En 



(86) 
(87) 



2. Spherical coordinates 

To study the problem in this coordinate system, we shall use the ex- 
pression of the partial Green's function (82) and use the following change of 
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variables 



{^^V) — *■ I v2rcos-, v2rsin- I . 
The partial Green's function (82) then can be written 

G^ir,ey,e';E) = f:(m+A, + i)I>^±^^p-ii^(cos^) 

da 



(88) 



m=0 



x^4m(cos^')t 



2' m 
2Mlo 



inVrPJo sin(a;cr) 



exp 



h 



{2a+i0)a 



X exp 



iMu! 



r + r') cot{uja) 



( 2Mu 



'2m+2Ai+l 



Y i^sin(a;cr) 



(89) 



We have used here the addition theorem (36) , the connection between hy- 
pergeometric functions and the Jacobi polynomials (see Ref.^^ , p. 1036, Eq. 
(8.962)) 



2Fi(/ + a + /3 + l,-/;l + a; 



l\r{a+l) M),. , . 



as well as the relation between the hypergeometric functions and the Legen- 
dre polynomials (see Ref.^^ , p. 1009, Eq. (8.771)), and eventually the hnk 
between P-"*(x) and P^{x) (see Ref.=^^ , p. 1008, Eq. (8.752.2)). 

In order to perform the integration on the a time variable, we shall use 
the formula (40). Consequently, the final expression of the Green's function 
in spherical coordinates will be 

^ ^, Y^'^/ X l,r(m+Ai + l)r(p+m+Ai+l) 



n=0 m=0 
1 



M 



lurr 



——r 1 W 
n 



r(2m+2Ai+2) 

f2Mu 



xP-^A, (cos 9)P-1\^ (cos e')$P-^^^) (0)$i^^-^^^) (00 , 

(91) 



where r > r' and P——^- 
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The Hille and Hardy formula (23) and an analytic proceeding consisting 
in using the Sommerfeld- Watson transformation^^ will help us to write the 
Green's function (91) in the form of a partial wave development consisting 
of two contributions of a discrete and a continuous part: 

^ ' ' ^ ] E + iO-EN 

n,m=Q [, 1=0 

+ / , fe2fc2 > , 92 

Jo E + iO-!^ J 

with-for bound states-the poles located around the values of E and the nor- 
malized wave functions respectively given by 

^^ = "2^' ^ = ^ + ^ + Ai + l, (93) 



*i,m,n(r, ^, 0) = 



m+Ai + |)/!r(m+Ai + l) 



a3(/+A2+|)^m!r(Z+2A2+l). 



2r 



\2-i 



X exp 



a(l + A2 + I) 



-2X2 



a(Z+A2+|)y 
2r 

a{l + A2 + I) 



x^^-i\,(cos^)$i±^-±'=^)(0). 



(94) 



For continuous states, the normalized wave functions and energy spectrum 
are respectively given by 



27r 

e 2ak 

r 



. ^ l,r(m + Ai + l) 



^ |r(| + A2 



ok 



r(2A2 + 1) 

x^^^,A.(-2^^r)p-i\^(cose)$(±^^'±^^)(0), (95) 



V. Three-dimensional minimally super-integrable po- 
tentials 

There are nine three-dimensional potentials which belong to the class 
of minimally super-integrable Smorodinsky-Winternitz potentials, that is to 
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say three-dimensional potentials characterized by the existence of four func- 
tionally independent integrals of motion. Among them are seven potentials 
which have S0(2, 1) as a dynamical group and thus their exact solution can 
be given via the (M-S) variant of the algebraic approach in different coordi- 
nate systems. As an example, we shall discuss the potential 



kfxs 



+ F 



X2 



(96) 



with ki a positive constant. The corresponding observables have the form 

+ F(tan0), 



11 

2M 



p2 



1 
4M 



' 2M ' 2M 

)-«ofg + 



T 

-"l 2M 

kf cos9+F(tan(; 

3-; 



sin" y ' 

^(i-i) (A;? + F(tan0)), 



(97) 



2M 



where lij = {Li, Pj} = UPj + PjLi with (i, j) = {xi, X2, x^) 

The Green's function for this potential can be explicitly evaluated in the 
parabolic and spherical coordinate systems. For F (|) = 7^ and 7 a real 
constant, this potential reduces to the ring-shaped potential proposed by 
Hartmann as a model for the ring-shaped molecules. It has been analyzed 
by many authors in the framework of path integrals^^"^'^ and through the 
algebraic approach^^"^^ .We can also notice the close link of the latter with 
the Coulomb potential plus the barrier created by the solenoid of Aharonov- 
Bohm^^ treated with the path integrals^°"^^ and via the algebraic technique^^ 
. To give the solution for the potential (96) via the (MS) variant of the so(2, 1) 
algebraic approach, we will use the Kustaanheimo-Stiefel transformation^^ 
{xi,X2, X3} {ujj e [1, 4]} corresponding to the surjection — > R^, which 
can be defined as: 



(A) 



with the constraint 



U2 
U3 
\U4 J 



(A) 



/ 


U3 


U4 


Ui 


U2 


\ 




—U4 


Us 


U2 


-Ml 






-Ui 


-U2 


U3 


U4 




\ 


-U2 


Ui 


—U4 


U3 


/ 



(98) 



dx4 — 2{—U2dui + Uidu2 — 1*4^^3 -|- ^3^1*4) = 0, 



(99) 
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allowing to define a fourth coordinate 

X4 — 2 J {—U2 ill +Ui U2 —114 U3 +Us ii^jds. 
Moreover, we can show that 



(100) 



/ 


d 
dxi 


\ 




/ 


8 


\ 




d 


















da. 






±± 
2r 8x4, 








J 


V 


) 




V 


5! 

8U4 



(101) 



The (KS) transformation allows us to write the Laplacian in i?^ in 
terms of the laplacian in i?^ as 



4r 



dxl ' 



4 02 

j=i J 



(102) 



where r — {xl + xl + xf) ^ — ul + ul + ul + ul- 

Using the Schwinger's integral representation, the Green's function rela- 
tive to the potential V4(~r^) is written 



G{r, r';E) - 



h 



2M 



6{ r — r' ). 

(103) 

It is possible to introduce an additional variable X4 by means of the well- 
known identity 



i: 



exp 



h \8Mr^dxly J 



5{x4)dX4 



and to show that expression (104) can be written as 

G(y,r' ;£;)=/ dxiGi'r' ,X4, r' ,0; E), 
J —00 

where 



(104) 



(105) 



G{l^,X4y,0;E) = dSexpl.-^-^[HT-E-iO]}6{V- r' )6{x4) 



(106) 
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with 



2M 



(107) 



Using (102), (98), the Jacobian of this transformation given by 

,3:2,3:3,3:4) _ iQj.2 ^Yie time transformation r — t-, the Green's function 

d(ui,U2,U3,U4) ' 4r' 

(106 ) can be put in the form 



1 

4r 



drexp 



h 



(4q;o + iO)T 



IT 

xexpi-- 



2M 



□2 + V{u) 



(108) 



where 



Viu) = 



2M Vitf + ii| ul + ulj 2M\ul + ul 



ul + ulJ \ U1U3 + U2U4 J 



(109) 



The evaluation of this expression is possible in two coordinate systems. 



3. Parabolic coordinates 

Going on to the double polar coordinates 



{ui, U2) {r], (pi) :ui^-^ cos 0i, 
{u3, U4) (^, 02) ■'^3^ cos 02, 



the Green's function ( 108 ) becomes 



G(T*, X4, r' , 0; £■) ~ 2~ I drexp 



U2 



^ sin 01, 
^ sin 02, 



-(2ao + iO)T 



exp 



(110) 



X^^^4pf^^(0.- 00^(0. -0'2), (111) 
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where 



with 



92 

+ 



Id Id 



1 



1 92 



2M \di^ drf ^d^ r]dr] dcpj rf 



(112) 



mi / 1 1 



2M 



^2 ny 2M VC^ 



- + - F(tan(0i-</.2)-^(e'+??'). (113) 



Here we have apphed the rescahng of 2t to r. 

At this point, we notice that the separation of variables is not possible. 
To achieve it, it is necessary to perform the integration on the variable X4 by 
using the Euler's angles 0i = 02 = (0 < < 27r, < a < 47r). 

Then, it is easy to see that dx^ = rda and that by integration over a we 
are led to 4:7r5^^o. It follows that the Green's function (104) can be written 



iT ^2 



GiT", r';E)= / drexp 
Jo 

where 

Ki^{(t),<p']T) = exp 

and 
with 



-{2aQ + iO)r 
a 



1 1 



K^,{<l>,<l>';r)K(^,rj,^',rj';r), 

(114) 



h 2M VC^ 77^ 



+ F(tan0) 



(115) 



1 



92 92 

+ 



= exp 



SiC-e)Siv-v'), (116) 



2M \d^'^ drj^ 



2M 



(117) 

In order to bring to a constant the mass appearing in the Hamiltonian 
contained the kernel to expression (115) let's perform the time transformation 

T — > (7 defined hy a — t (^-^ + . Then, it follows that 



= exp 



T2M 



92 



+ F(tan0) 



5(0 - 00 
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dEx^ exp 



(118) 



with Ex^ = '■^2i4'- Here, we have assumed that the propagator associated with 
the potential F{tan(f)) is known. 

Let's now insert (118) and (116) in the expression (114), we obtain 

1 



G'( r , r';E) 



where 



xK((,(';T)K{r,,!ri';T), 



dr exp 



h 



{2ao + iO)T 



(119) 



K{u,u';t) — exp <^ — — 



XI T kf 



2 ^ U2 1 
4 



S(u - u') 



exp 



IT 

T 



2M \dv? 

[T^{u) + 2hWT^{u)\ )5{u-u')-{u = i,r]), (120) 



with u = \/-^. 



Following formula (13), the kernels (120) are written 



K{u, u'; t) 



MU! 



exp 



iMuj 

~2h^ 



u -\-u ) cot(a;T) 



Mojuu' 
ihsmiujr) 



(121) 

with = ^ kl. 

Inserting (121) into (119) and by a procedure similar to that which has 
led to result (83), we obtain 

^ ' ' ^ J M E + iQ-Ej, 

+i, — 

with the normalized wave functions and the energy spectrum for 



, (122) 
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-the bound states: 



*ni,n2,A^(C,^,0) 



ni\n2\ 



a^m r(ni + A_ + l)r(n2 + A+ + 1) 



2/^X2 

aN 



X*A,(0), 



\ rA- / ^ \ rA+ 



(123) 



^iv = -^^, iV = ni+n2 + -(A_ + A+) + l, (124) 
-the continuous states: 

r (1(1 + A+) + 1 (i + ^)) r (i(i + A-) + i(i - 



*P,.,A,(e,?7,0) 



27rr(l + A+)r(l + A_) 

g 2ap 

y^^?7 -2F(.a+'^^'^ ~2iia~^)'^ 

X*A, (</>)■ 



2^2 



hp 

2M 



(125) 
(126) 



4. Spherical coordinates 

With the help of the change of variables defined by 

(C,??) (V2rcos^,y2rsin^), 



(127) 



and by applying the Bateman 's expansion formula (26), expression (119) is 
rewritten 



/oo 
^£^A,*A,(0)*l,(0')E*n 
n=0 

/sin 9 sin ^' Jo 



^(A+,A_) 



Mu 

dr- 



ihsm(ujT) 
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X exp 



-{2ao + i0)r 



exp 



(r + r') cot(a;T) 



/ 2MuVrr' 

2n+A++A_+l I -t . / T 



(128) 



Performing the integration over the time variable r with the help of for- 
mula (30), we obtain the following final form 



G{l',r';E) = / ci£;,,*,,(0)*^^(0') E T^^^*" ^' 
J ^ V sm sm B' 



A_) / I 



1 r(p + n + i(A+ + A_) + 1; 



2 7 2ia;rr' r(2n + A+ + A_ + 2) 

X^-p,n+i(A++A_+l) I 



where r > r' and P = — ^ ■ 



Following the calculation procedure in subsection (IV. 2), one shows that 

(130) 



n=0 I, i=0 







where the normalized wave functions and the energy spectrum are given by 

Al-i 



*/,n,A.(^) 



a (; + Ai + I) 
X exp 



2l\ 



ar(; + 2Ai + l) 



2r 



a (/ + Ai + I) 
x$i-^-'-n^)*A,(0), 



2Ai 



a (Z + Ai + I) 
2r A / 2 



(/ + Ai + i)y V sin^ 
(131) 
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Ma"^ 1 1 

^Ar = -2^, N^l + X, + - and Ai = n + -(A+ + A_ + 1), (132) 

for the bound states and 

x*A,(</>), (133) 



^fe = (134) 
fe 2M' ^ ^ 



for the continuous states. 



VI. Conclusion 

In this paper, we have analyzed through the Milshtein and Strakhovenko 
variant of the so(2, 1) Lie algebra a set of potentials belonging to three dif- 
ferent classes of Smorodinsky-Winternitz potentials. The use of the second 
order differential operators of this algebra allows to write the Hamiltonian 
of these physical systems in form of a linear combination of the latter. Us- 
ing the Schwinger's integral representation and with the help of two Baker- 
Campbell-Hausdorff formulas allowing the separation of the Tj operators and 
thus simplifying their action on a Laplace transform of a well chosen Dirac 
distribution , we have shown that we can construct the Green's functions in 
compact form in different coordinate systems. This method can be compared 
to the approach of the Schrodinger equation and to the Feynman's path in- 
tegral technique. It gives a local view of the problem under consideration 
like the Schrodinger approach, but its advantage is in the computation of 
the explicit and compact form of the Green's function from which the energy 
spectrum and the suitably normalized wave functions are simultaneously ex- 
tracted for the bound states and for the continuous states if they exist at one 
and the same time. 

The advantage of the path integral approach in comparison with the 
algebraic method is that it provides a global view of the dynamics of the 
physical system, but a problem of singularity is often to be found at the origin 
of coordinates, and it requires a regularization which is rather complicated to 
perform. Prom this point of view, we can assert that this algebraic method 
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has the advantage of not presenting this problem owing to the fact that it is 
locaL 

The method of Milshtein and Strakhovenko could become a powerful 
alternative approach to the path integral technique if we manage to extend 
its use to the treatment of Poschl- Teller potential class. 



31 



References 

^H. Kleinert, Fortschr. Phys. 6, 1 (1968); Group Dynamics of the Hydrogen 
( Lectures presented at thel967 Boulder Summer School, published in 
Lectures in Theoretical Physics, Vol.XB, A. O. Barut and W. E. Brittin, 
eds. Gordon and Breach, New York, 1968 ). 

^G. A. Natanzon, Teor. Mat. Fiz. 38, 146 (1979). 

^M. J. Englefield, Group Theory and the Coulomb Problem ( Wiley, New- 
York, 1972 ). 

4J. N. Huffaker and P. H. Dwivedi, J. Math. Phys. 16, 862 (1975). 

5 J. Cizek and J. Paldus, Int. J. Quantum Chem. 12, 875 (1977). 

^M. Berrondo and A. Palma, J. Phys. A: Math Gen. 13, 7 (1980). 

^J. Wu and Y. Alhassid, J. Math. Phys. 31, 557 (1990). 

8p. Cordero and S. Salamo, J. Phys. A: Math. Gen. 24, 5299 (1991); Foun. 
Phys. 23, 675 (1993). 

'^I. L. Cooper, J. Phys. A: Math. Gen. 26, 1601 (1993). 

^°C. C. Gerry, Phys. Lett. A 118, 445 (1986). 

i^H. Hartmann, Theor. Chem. Acta 24, 201 (1972). 

i^M. V. Carpio-Bernido and C. C. Bernido, Phys. Lett. A 137, 1 (1989). 

i^C. Quesne, J. Phys. A: Math. Gen. 21, 3039 (1988). 

^''A. I. Milshtein and V. M. Strakhovenko, Phys. Lett. A 90, 447 (1982). 

^^L. Chetouani, L. Guechi and T.F. Hammann, Helv. Phys. Acta 65, 
1069 (1992); Czech. J. Phys. 42, 913 (1992). 

i^A. N. Vaidya and H. Boschi-Filho, J. Math. Phys. 31, 1951 (1990). 

^^H. Boschi-Filho and A. N. Vaidya, Ann. Phys., (NY) 212, 1 (1991). 



32 



i^H. Boschi-Filho, M. de Souza and A. N. Vaidya, J. Phys. A: Math. Gen. 

24, 4981 (1991). 

^^L. Chetouani, L. Guechi, T. F. Hammann and A. Lecheheb, J. Math. Phys. 
36, 1016 (1995). 

Chetouani, L. Guechi and T. F. Hammann, J. Math. Phys. 33, 3410 

(1992) . 

2iH. Boschi-Filho and A. N. Vaidya, Phys. Lett. A 149, 336 (1990). 

^^L. Chetouani, L. Guechi, A. Lecheheb and T. F. Hammann, Helv. Phys. 
Acta 68, 158 (1995). 

^^A. A. Makarov, J. A. Smorodinsky, Kh. Valiev and P. Winternitz, 
Nuovo Cimento A 52, 1061 (1967). 

24m. Kibler and P. Winternitz, J. Phys. A: Math. Gen.20, 4097 (1987); Phys. 
Lett. A 147, 338 (1990). 

^^N. W. Evans, Phys. Lett. A 147, 483 (1990). 

26n. W. Evans, J. Math. Phys. 31, 600 (1990); ibid 32, 3369 (1991). 

2^D. Bonatos, C. Daskaloyannis and K. Kokkotas, Phys. Rev. A 48, R 3407 

(1993) ; J. Hietarinta, Phys. Rep. 147, 87 (1987). 

W. Evans, Phys. Rev. A 41, 5666 (1990). 

2^C. Grosche, G. S. Pogosyan and A. N. Sissakian, Fortschr. Phys. 43, 
453 (1995). 

^°H. Klcincrt, Path Integrals in Quantum Mechanics, Statistics and Polymer 
Physics ( World Scientific, Singapore, 1990 ). 

^^R. Gilmore, Lie Groups, Lie Algebras and Some of Their Applications 
( Wiley, New York, 1974 ). 

^2 J. Schwinger, Phys. Rev. 82, 664 (1951). 

^^B. G. Wybourne, Classical Groups for Physicists (Wiley, New York 1974). 

33 



H. Duru and H. Kleinert, Fortschr. Phys. 30, 401 (1982). 
Inomata, Phys. Lett. A 87, 387 (1982). 
^^T. Levita Civita, Opere Matematiche, V., 2 (1956). 

^^I. S. Gradshtein and I. M. Ryzhik, Table of Integrals, Series and Products 
( Academic, New York, 1965 ) p. 1038, Eq. (8.976). 

^^A. Erdclyi, W. Magus, F. Oberhettinger and F. G. Tricomi, Higher Tran- 
scendental Functions ( Mc Graw-Hill, New York, 1953 ) Vol. 2. 

^^P. M. Morse and H. Feshbach, Methods of Theoretical Physics ( Mc Graw- 
Hill, New York, 1953 ) Vol. 1, p. 781. 

^°H. Buchholz, The Confluent Hypergeometric Function, Springer Tracts in 
Natural Philosophy ( Springer- Verlag, Berlin, 1969 ) vol.15. 

^^I. H. Duru and H. Kleinert, Phys. Lett. B 84, 185 (1979). 

^2r. Ho and A. Inomata, Phys. Rev. Lett. 48, 231 (1982). 

^^K. Pak and L Sokmen, Phys. Lett. A 100, 327 (1984). 

^^A. Inomata, Phys. Lett. A 101, 253 (1984). 

^^F. Steiner, Phys. Lett. A 106, 363 (1984). 

^^L. Chetouani and T. F. Hammann, J. Math. Phys. 27, 2944 (1986). 
^^L. Chetouani and T. F. Hammann, Nuovo Cimento B 98, 1 (1987). 
^^H. Kleinert, Phys. Lett. A 120, 361 (1987). 

^^D. P. L. Castrigiano and F. Stark, J. Math. Phys. 30, 2785 (1989). 

5°S. N. Storchak, Theor. Math. Phys. 82, 32 (1990). 

^^C. Grosche, Fortschr. Phys. 40, 695 (1992) . 

^^M. V. Carpio-Bernido, J. Math. Phys. 32, 1799 (1991). 

s^M. Kibler and C. Campigotto, Int. J. Quantum Chem. 45, 209 (1993). 

34 



^''M. Kiblcr, L. G. Mardoyan and G. S. Pogosyan, Int. J. Quantum Ghem. 
52, 1301 (1994). 

^^W. Thirring, Classical Field Theory ( Springer- Verlag, New York, 1986 ) 
p. 140. 

^^M. V. Garpio-Bernido, G. G. Bcrnido and A. Inomata, in Third Interna- 
tional Conference on Path Integrals from meV to MeV, 1989, p. 442, Eds. 
V. Sa-Yakanit et al ( World Scientific, Singapore, 1989 ). 

^^M. V. Garpio-Bernido and A. Inomata, in Bielefeld Encounters in Physics 
and Mathematics VII, Path Integrals from meV to MeV, 1985, p. 261, Eds. 
M. G. Gutzwiller, A. Inomata, J. R. Klauder and L. Streit ( World Sci- 
entific, Singapore, 1986 ). 

^^I. Sokmen, Phys. Lett. A 115, 249 (1986). 

^^L. Ghetouani, L. Guechi and T. F. Hammann, Phys. Lett. A 125, 277 
(1987). 

™M. V. Garpio-Bernido, J. Phys. A: Math. Gen. 24, 3013 (1991). 

^^M. Kibler and G. Gampigotto, Phys. Lett. A 181, 1 (1993). 

^^M. Kibler, G. H. Lamot and P. Winternitz, Int. J. Quantum Ghem. 43, 
625 (1992). 

6^F. Galogero, J. Math. Phys. 10, 2191 (1969). 

6^M. Kibler and T. Negadi, Int. J. Quantum Ghem. 26, 405 (1984). 

6^A. Guha and S. Mukherjee, J. Math. Phys. 28, 840 (1987). 

^^Ya. A. Granovsky, A. S. Zhedanov and I. M. Lutzenko, J. Phys. A: Math. 
Gen. 24, 3887 (1991). 

6^A. S. Zhedanov, J. Phys. A: Math. Gen. 26, 4633 (1993). 

^^M. Kibler, L. G. Mardoyan and G. S. Pogosyan, Int. J. Quantum Ghem. 
52, 1301 (1994). 



35 



'Y. Aharonov and D. Bohm, Phys. Rev. 115, 485 (1959). 

'L. Chetouani, L. Guechi and T. F. Hammann, J. Math. Phys. 30, 655 
(1989). 

I. Sokmen, Phys. Lett. A 132, 65 (1988). 

'Lc Van Hoang and N. T. Giang, J. Phys. A : Math. Gen. 26, 3333 (1993). 

'M. Kibler and T. Negadi, Phys. Lett. A 124, 42 (1987). 

P. Kustaanheimo and E. Stiefel, J. Reine Angew. Math. 218, 204 (1965). 



36 



